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Abstract 

The relativistic perturbative approach based on the Bethe-Salpeter (BS) equation 
is developed for the study of the characteristics of the hadronic 7r + 7r~ atom. The 
general expression for the atom lifetime is derived. Lowest-order corrections to the 
relativistic Deser-type formula for the atom lifetime are evaluated within the Chiral 
Perturbation Theory. The lifetime of the 7r + 7r~ atom in the two-loop order of the 
Chiral Perturbation Theory is predicted to be T\ = (3.03 ± 0.10) x 10~ 15 s. 
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The study of the pion-pion scattering process enables one to gain a deeper insight in 
the nature of strong interactions. According to common belief, low-energy pion dynamics 
is described in the Chiral Perturbation Theory (ChPT) ^ which predicts the 5-wave tttt 
scattering lengths a[j = 0.217m" 1 and = —0.041m" 1 with an 5% accuracy 0. The 
calculations within the Generalized ChPT, involving more parameters, lead to a most 
likely values of 7T7T scattering lengths || which significantly differ from the predictions of 
the standard ChPT (e.g. a[j = O^Tm^ 1 ). Despite this difference, both results for the 
scattering length a{] are compatible with the experimental value a|] = (0.26±0.05)m^ 1 [||] 
due to a large experimental error. Consequently, a precise measurement of 7T7T scattering 
lengths will be an excellent test of the ChPT. 

The existing experimental information on the S-wave ttit scattering lengths is extracted 
from the data on the process ttN — > tttt N and on K e ^ decays [|J . The DIRAC experiment 
at CERN is aimed at the precise measurement of the lifetime of the 7r + 7r~ atom This 
will allow to determine the quantity a{] — clq with a high accuracy and thus will provide a 
decisive probe of the predictions of ChPT. The experiments on the measurement of the 
hadronic atom observables are held also at PSI, KEK, Frascati, Uppsala. 

The experimental study of 7r + 7r~ atom characteristics provides a direct information 
about the strong tttt amplitudes near threshold since the average momenta (~ am,) of 
constituents in the atom is much less than either of two typical scales characterizing 
strong interactions in chiral theories (~ m n and 4irF n ~ 1 GeV). The huge difference 
between the bound-state and strong scales results in a clear-cut factorization of strong and 
electromagnetic interactions in the observables of the pionium: the Coulomb interaction 
between pions is mainly responsible for the atom formation, and strong interactions lead 
to its decay. In the lowest-order formula for the ground-state atom lifetime |J these two 
distinct phenomena enter separately through the value of the nonrelativistic Coulomb 
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ground-state wave function at the origin ^i(O) and the difference of the S-wave scattering 
lengths aj] — Oq, respectively. 
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where, using the values of the scattering lengths from the two-loop calculations in standard 
ChPT a ° -a 2 = 0.258m; 1 , we obtain n = 3.20 x 10" 15 s. 

In the view of the forthcoming DIRAC high-precision experiment which is aimed at the 
extraction of the difference a{j — from the pionium lifetime measurement, it is important 
to estimate the precision of Eq. (H) itself and calculate all possible corrections to it which 
were neglected in the original derivation. In particular, one can readily observe from 
Eq. (0) that the ad hoc change of the "reference mass" from the charged to the neutral 
pion mass in the definition of strong scattering lengths leads to a ~ 12% variation in the 
lifetime. The leading electromagnetic corrections of order 0(a) will also be important in 
the analysis of the high-quality data from DIRAC experiment. 

The earlier efforts of the evaluation of the corrections to the pionium lifetime in the 
field-theoretical framework |7j] were focused on the evaluation of the quantities entering 
into Eq. ([I]). These investigations resulted in the understanding of the fact that the 
consistent approach to the calculation of bound-state characteristics can not be confined 
solely to the corrections to the scattering lengths and atomic wave function. Instead, 
to this end a systematic field-theoretical framework is needed to consistently take into 
account the pion strong dynamics described by the effective ChPT Lagrangian in order to 
provide a reliable and unambiguous test of the predictions of ChPT. In the several recent 
papers an attempt was made to carry out a comprehensive field theoretical analysis of 
the pionium problem. Namely, in || on the basis of the 3-dimensional constraint theory 
bound-state equations the mass shift, radiative and second-order strong corrections have 
been evaluated. In our previous papers || a systematic perturbative approach based on 
the BS equation is constructed, and a complete set of the lowest-order corrections to the 
pionium lifetime is evaluated. The vacuum polarization effect on the pionium observables 
is evaluated in within the nonrelativistic QED. 

In a nonrelativistic potential approach to the problem (Tl| the expression for the life- 



time, including the corrections, is written in terms of physical scattering lengths in charged 
and neutral channels in the presence of Coulomb interactions. This is a counterpart of the 
lowest-order Deser-type factorization: the bound-state observables are expressed in terms 
of the quantities characterizing the scattering process, and these quantities differ from the 
purely strong ones by the mass shift and electromagnetic corrections. In order to evaluate 
these corrections, one needs the titt potential to be given explicitly. Assuming local and 
energy-independent strong potentials in the isospin symmetry limit, one reproduces the 
7T7T scattering phase shifts calculated from ChPT. At the next step the Coulomb potential 
is added and the mass shift effect is included via assigning the physical masses to pions, 
the strong potential remaining unchanged. 

Below we give the basic ideas and assumptions of our approach to the pionium ob- 
servables, based on the BS equation. 

The pions within our approach are described by the elementary fields in the La- 
grangian. This allows one to write down the bound-state BS equation for the metastable 
7T + 7r~ atom. The necessary link to the strong amplitudes from ChPT where the pions are 
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described by the pseudoscalar quark densities, proceeds via the Deser-type factorization 
in analogy with the potential theory. Namely, we express all corrections to the Deser 
formula in terms of the on-mass-shell pion amplitudes which are further identified with 
the amplitudes from ChPT. 

In our calculations we always adopt the so-called "local" approximation which consists 
is suppressing the relative momentum dependence of the strong amplitudes. The origin 
of this approximation can be traced back to the huge difference between the bound-state 
and strong momentum scales. 

The pionium in our approach is described by the exact bound-state BS equation with 
the full kernel V(P) which, apart from all two-particle irreducible diagrams includes the 
self-energy corrections in two outgoing charged pion legs ||. The square of the cm. 
momentum P takes the value P 2 = M 2 = M 2 — iMT where M denotes the "mass" of an 
atom, and V stands for the decay width. 

In order to perform the perturbative expansion of the bound-state observables we 
single out from the full kernel V the instantaneous Coulomb part Vc which is responsible 
for the formation of the bound state composed of ir + and ir~ . The "remainder" of the 
potential which is treated perturbatively is denoted by V = V — Vc- 

The "unperturbed" part of the full kernel Vc is chosen in the form guaranteeing the 
exact solvability of the corresponding BS equation. The solution of this equation yields the 
relativistic Coulomb wave function ipc with the corresponding eigenvalue P* 2 = M* 2 = 
m 2 (A — a 2 ) where m n is the charged pion mass. Then, the solution of the full BS equation 
is expressed in terms of the unperturbed solution ||. Substituting this solution into the 
complete BS equation we arrive at the relation 

< lM[l + (Go\P) ~ G \n - V'{P))G R Qr\G^\P) - G Q \n - V'(P)Mc >= 

which serves as the basic equation for performing the perturbative expansion of the bound- 
state observables. Here Go denotes the free Green's function of the tt + tt~ pair, Gr stands 
for the pole-subtracted part of the exact relativistic Coulomb Green's function and Q is the 
projection operator on the subspace orthogonal to the unperturbed ground-state BS wave 
function ||. The explicit form of the operator GrQ is given and the perturbation kernel V 
is known up to any given order of loop expansion. The only unknown quantities entering 
into this expression are the real and imaginary parts of P 2 . Consequently, expanding this 
expression in powers of V, we arrive at the recursive relations which define the mass and 
width of the bound state in any given order. 

In order to provide the necessary Deser-type factorization in the obtained results, one 
needs a further classification of diagrams entering into the definition of V. We consider 
the following subsets of diagrams: 

1. A purely strong part, which is isotopically invariant. 

2. The part which is responsible for the m^i — m^o electromagnetic mass difference 

3. Remaining electromagnetic effects, including the exchanges of virtual photons. 

The sum of the potentials corresponding to the parts 1 and 2 is denoted by Vi 2 = 
V\ + Vi. The T-matrix corresponding to summation of the potential Vu in all orders is 
given by Ti 2 (P) = Vu(P) + Vi2(P)Gq(P)Ti2(P) . The rest of the potential is referred to 
as V3 = V — V12 and is treated perturbatively. Further, the quantity GrQ in our basic 
relation can be written as GrQ = G (P*) + 5G where 5G corresponds to the exchange 
of the ladder of Coulomb photons and is also considered as a perturbation. Then in the 
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first order in V3 and 5G the basic relation takes the form 



= -2iM*5M- < ipc\Tn\^c > + 

+ 5M< ^c\{G^)'G Q T 12 \^ c > +\{5M) 2 < VcKQ 1 )" + 2T 12 G (G 1 )"\ij c > + 

+ < VcK^Gq 1 )' - T 12 )5GT 12 \ip c > - < ip c \(l + T 12 G )V 3 (1 + G T 12 )\^ C > (2) 

where Go = Gq(M*) and the prime stands for the differentiation with respect to M*. 

In the "local" approximation T\ 2 does not depend on the relative momenta and thus 
can be taken out of the matrix elements. Neglecting first the corrections (only the first 
line of Eq. (0) contributes), we arrive at the relativistic counterpart of Eq. (p]) 

-1/2 EM } = 2M^ |Vl/l(0)l { Im } W-^r- «; 0, g0) ) (3) 

where T(s;p,q) = 167rT 12 denotes the 5-wave 7T7t strong scattering amplitude which 
includes the effect of m„-± — m n o mass difference. g*o is the relative momentum of the 7r°7r° 
pair at the threshold s = 4m^, with the magnitude given by the relation m\ = m 2 ^ + q Q 2 . 

In order to calculate the corrections to the relativistic Deser formula, one has to 
evaluate the integrals entering into the Eq. (Q). 

1. The correction due to the shift of the bound-state pole by strong interactions is 
determined by two terms in the second line of Eq. (§). Calculating these integrals explicitly 
and using Eqs. (|3]) we obtain 

5 S = -9A£ (1) /(8£i) = -5.47 x KT 3 m w (2a° + a 2 ) (4) 

where E\ denotes the unperturbed ground-state binding energy. 

2. The correction due to the relativistic modification of the Coulomb wave function is 
given by the value of the wave function ipc & t the origin which for the particular choice 
of the kernel used in the calculations is related to its nonrelativistic counterpart by 

V>c(0) =*i(0) (l-0.381a + ---) (5) 

and the correction in the decay width which is proportional to |-?/>c(0)| 2 is twice as large. 

3. The correction due to the exchange of Coulomb photon ladders is given by the first 
term in the third line of Eq. (§) containing 5G. Using the known expression for the 
relativistic Coulomb Green's function and evaluating this integral explicitly, we obtain 

5 C = (1/2 + 2.694 - lna) AE (1) /E 1 = 3.95 x 1(T 2 m n (2a° + a 2 ) (6) 

4. In order to obtain the mass shift and radiative corrections, in the perturbation kernel 
V3 one has to take into account the residual photon exchange diagram and the self-energy 
corrections in the outgoing charged pion legs. The net result of this effect consists in the 
cancellation of the correction due to the relativistic modification of the Coulomb wave 
function and in replacing the quantity 7^ r + 7r -_ >7r o 7r o which appears in the imaginary part 
of the 7T7T scattering amplitude in Eq. (||), by the Coulomb pole removed full amplitude 
for the process it + 7t~ — > 7r°7r°. Using the expression of this amplitude |12[], and isolating 



the purely strong isotopically symmetric amplitude with a common mass equal to the 
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charged pion mass, we arrive at the following expressions for the mass shift and radiative 
corrections 



>M 



^( 1 + fi3fe(«^-T^*-«')) < 7 » 



S em = a/(12rr) (-30 + 3/Cf - ICf) (8) 

where /Cf and U denote the low-energy constants of ChPT |lj, |12| . 

5. Taking into account the diagram in the kernel V3 corresponding to the vacuum polar- 
ization effect and evaluating numerically the resulting integral, we arrive at the following 
result 

5 vac = 3a 2 m n /{16 m e ) x 0.6865 (9) 

which completely agrees with the zero-Coulomb piece of the result given in Ref. [|Kj . 

6. The correction to the pionium decay width due to the finite size effect is caused by the 
modification of the instantaneous Coulomb interaction by the pion loop in the two pion- 
photon vertex in the kernel V3. For simplicity, we approximate this vertex by a monopole 
parameterization. After calculating the corresponding matrix element explicitly, we obtain 

5 F = laml < r 2 > v /(3tt) In {m\ < r 2 > v /24) (10) 

where < r 2 > v denotes the square charge radius of the pion. 

The relativistic Deser formula with all above corrections takes the form 

1 167r /2Am 7r \ 1 / 2 / i Am^ 1 / 2 , n 2 . 2 . _ , n ..„ .., . . . . . . . 

T i = — \~^^) ( a o~ a o) l*i(0)| 2 (l+Ss+5c+SM+S em +5 V ac+5 F ) 

where a{] and denote the tttt scattering lengths in the isospin-symmetric case, with the 
charged pion mass taken to be the common mass of the pion isotriplet. 

In our calculations for the constants U we take the numerical values from Ref. [p]] 
h = -2.3 ± 2.7, l 2 = 6.0 ± 1.3, l 3 = 2.9 ± 2.4, f 4 = 4.3 ± 0.9. Also, we use the values 



/Cf u = 1.8 ± 0.9, t^f-Kf 1 = 0.5 ± 2.2 from Ref. |J. Also, we adopt the following 



value of the e.m. charge radius of pion < r 2 >y= 0.439 Fm [|l]. With the use of the 
above values of the parameters the total correction to the pionium decay width turns out 
to be5 tot = (6.1 ±3.1)% and 

ri = (3.03 ±0.10) x 10~ 15 s 

The largest correction in the decay width is caused by the total effect of the mass shift and 
electromagnetic radiative corrections. The recent field-theoretical calculations carried in 
the framework of the 3D relativistic equations have come to a general agreement with our 
results for all individual effects. However, the sign of the mass shift effect obtained in the 
nonrelativistic scattering theory approach [[□]] turns out to be opposite as compared to 



our result, and is of the same order of magnitude. Owing to the derivative character of 
pion couplings in the chiral Lagrangian, a possible reason for this discrepancy might be 
an explicit energy independence of strong potentials used in these calculations. 

To summarize, we have evaluated a complete set of the lowest-order corrections to the 
pionium decay width in the "local" approximation for the strong tctc amplitude. For a 
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full understanding of the problem, however, the reason for the difference in sign in the 
mass shift effect in the field-theoretical and potential theories should be investigated in 
detail. Also, it is very important to have a reliable quantitative estimate of the accuracy 
of "local" approximation including a proper treatment of arising new UV divergences. 
Note that recent attempts of taking into account the momentum dependence of strong 
amplitudes have so far led to a controversial results for the correction to the decay width 
~ 8% |L3|], ~ 2% |L4| (large values of the pionium lifetime obtained in the latter paper 



are merely due to the use of the tree- level scattering lengths in the Deser formula). 
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